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NOMENCLATURE 


Cross  sectional  area  of  stiffeners 
Stiffener  spacing 

Young's  modulus  of  shell  and  stif¬ 
fener,  assumed  same 

Stiffener  eccentricity  (dimensionless 

Dimensionless  stress  function 

Shear  modulus  of  skin  and  stiffener 

Shell  thickness 

Moment  of  inertia  of  stiffener  about 
shell  middle  surface 

Torsional  constant  of  stiffener 

Cylinder  length 

Wave  numbers  in  axial  and  circum¬ 
ferential  directions 

m  ,  m  )  Stress  couples  (dimensionless) 
xy  yx 

Stress  resultants  (dimensionless) 

Critical  load  (dimensionless) 

Radius  of  shell 

Middle  surface  displacements 
(dimensionless) 

Axial  and  circumferential  coordinates 
(dimensionless) 

Batdorf  parameter  =  L2/RH 

Axial  stiffness  parameter  = 


111 


ma 


s  ’ 


V 

s,r 

comma 


Flexural  stiffness  parameter  = 


Torsional  stiffness  parameter  = 


Poisson's  ratio  of  shell  and  stiffeners 

Denote  stringers  (axial  stiffeners) 
and  rings 

Denotes  partial  differentiation 


INTRODUCTION 


Research  activity  in  the  compressive  general  instability  of  stif¬ 
fened  cylindrical  shells,  in  spite  of  some  early  work  in  the  1930 's 
(see  for  example,  References  1  and  2),  did  not  receive  sustained  effort 
until  the  advent  of  launch  and  space  vehicle  designs  of  the  early  60' s. 

In  this  latter  period,  extensive  test  results  were  obtained  on  stiffened 
cylinders  which  showed  that,  by  and  large,  these  agreed  with  the  predic¬ 
tions  of  a  Donnell-type  linear  theory  which  included  the  stiffeners  in 
the  expressions  for  the  resistance  of  the  shell. 

Of  the  various  formulations  used  in  this  period,  we  are  concerned 
here  with  the  linear  theory  of  Prof.  Singer  and  his  coworkers  at  Technion,3 4 
which  was  also  proposed  by  Block  et  a 2l1*  in  a  slightly  more  generalized  form. 
The  main  contribution  of  this  theory  is  that  it  accounts  for  the  coupling 
between  the  bending  stress  resultants  and  the  middle  surface  stretching 
on  the  one  hand  and  that  between  the  membrane  stress  resultants  and  the 
changes  in  the  middle  surface  curvature  on  the  other.  This  coupling, 
produced  by  the  eccentricity  of  the  stiffener  center  line  location  with 
respect  to  the  shell  middle  surface,  has  the  net  effect  of  involving  odd 
powers  of  the  eccentricity  parameter  in  the  final  expression  for  the 
buckling  stress.  Not  surprisingly  then,  the  sign  of  the  eccentricity 
with  respect  to  the  middle  surface  of  the  shell  (i.e.,  due  to  the  stiffener 
being  outside  or  inside  with  respect  to  the  shell)  has  a  profound  influence 

1.  DSCUOU,  D.  Die  Druckfestigkeit  Versteifter  Zylindnscher  Schalen.  Luftfahrtforschung,  v.  ll.no.  8.  1935, 
p.  233-234. 

2.  TAYLOR,  J.  L.  The  Stability  of  Monocoque  in  Compression.  British  R  and  M,  no.  1679,  June  1935. 

3.  BARL'CH,  M.,  and  SINGER,  J.  Effect  of  Eccentricity  of  Stiffeners  on  the  General  Instability  of  Stiffened  Cylin¬ 
drical  Shells  Under  Hydrostatic  Pressure.  J.  Mech.  Eng.  Science,  v.  5,  no.  1 ,  1963,  p.  23-27 . 

4.  BLOCK,  D.  L.,  CARD,  M.  F.,  and  M1KULAS,  M.  M.  Buckling  of  Eccentrically  Stiffened  Orthotropic  Cylinders. 
NASA  TND-2960,  August  1965. 


on  the  buckling  load.  The  change  in  the  buckling  load  is  accentuated 
the  most  in  the  case  of  axially  stiffened  cylinders  under  axial  compression. 
This  eccentricity  phenomonon  is  discussed  extensively  by  Singer  et  al  in 
References  5  and  6. 

However,  since  their  main  preoccupation  has  been  with  this  eccentricity 
effect,  certain  aspects  of  the  problem,  such  as  the  characteristic  buckling 
modes,  are  not  adequately  treated  in  their  paper.5 6  For  instance,  in  the 
case  of  an  axially  stiffened  cylindei  the  axisymmetric  modes  are  dismissed 
with  the  statement  that  they  are  possible  for  certain  combinations  of 
the  geometrical  parameters.  Now,  it  is  a  characteristic  of  the  stiffened 
cylinder  buckling  that  unlike  the  unstiffened  c;  lenders  (where  the  linear 
solution  is  obtained  without  reference  to  the  buckling  modes)  the  solu¬ 
tions  are  governed  by  the  characteristic  modes:  axisymmetric ,  general 
asymmetric  and  for  axially  stiffened  cylinders,  an  asymmetric  mode  with 
a  single  half-wave  along  the  length  of  the  cylinder.7  In  fact,  for  some 
of  the  results  presented  in  Reference  6  pertinent  to  the  axial  case,  the 
buckling  mode  is  indeed  the  last  mentioned  type  with  a  single  axial 
half-wave,  though  nowhere  is  it  mentioned  in  the  paper. 

The  ab  initio  selection  of  the  governing  buckling  mode  for  a  given 
stiffener  geometry  is,  perhaps,  not  very  critical  in  obtainint  an  answer 
for  the  linear  problem;  however  it  becomes  very  important  in  a  related 
nonlinear  problem  connected  with  the  dynamic  stability  of  imperfect 

5.  SINGER,  J.,  BARUCH,  M.,  and  HARARI,  O.  Further  Remarks  on  the  Effect  of  Eccentricity  of  Stiffeners  on  the 
General  Instability  of  Stiffened  Cylindrical  Shells  TAE  Report  No.  42.  Technion  ■  Israel  Inst  of  Tech.,  Haifa. 
Israel,  Aug.  1965. 

6.  SINGER,  J.,  BARUCH,  M.,  and  HARARI,  O.  On  the  Stability  of  Eccentrically  Stifjened  Cylindrical  Shells  Under 
Axial  Compression.  Int  J.  Solids  Structures,  v.  3,  1967,  p.  445-470. 

7.  LAKSHMIKANTHAM,  C.,  GERARD,  G.,  and  MILLIGAN  R.  Linear  Orthotropic  Stability  of  Shallow-Stiffened 
Cylinders:  An  Assessment.  The  Aeronautical  Quarterly,  v.  2 1 ,  May  1970,  p  133-144. 
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cylinders.  Roth  and  Klosner8  have  demonstrated  that  in  the  case  of 
dynamic  stability  of  unstiffened  imperfect  cylinders,  the  static  linear 
(classical)  value  forms  the  limiting  solution  for  zero  imperfection. 

In  a  similar  study  with  stiffened  cylinders  the  present  authors  have 
found9'10  that  the  mode  of  buckling  remains  the  same  in  the  linear  (pel feet) 
and  nonlinear  (imperfect)  ranges,  and  hence  for  an  efficient  determination 
of  the  critical  dynamic  buckling  load  (that  is,  from  a  computer-time 
point  of  view),  a  knowledge  of  the  linear  static  modes  is  extremely 
important. 

Having  thus  established  a  rationale  for  reexamining  the  linear 
stability  problem  of  an  eccentrical ly  stiffened  cylindrical  shell,  we 
propose  to  show  that  from  among  the  numerous  geometrical  parameters 
entering  the  problem,  one  can  rationally  select  a  few  which  delineate  the 
regions  of  prevailing  buckling  modes. 

GOVERNING  EQUATIONS 

Although  the  buckling  equations  are  similar  to  those  of  Singer  et 
al,6  our  derivation  of  them  is  different,  being  patterned  on  the  non¬ 
linear  formulations.  Also  our  final  equations  are  cast  into  a  form 
which  lends  itself  to  an  easier  examination  of  the  relative  importance 
of  the  various  stiffening  parameters. 


8.  ROTH.  R.  S.,  and  KLOSNER,  F.  Nonlinear  Response  of  Cylindrical  Shells  Subjected  to  Dynamic  Axial  Loads. 
A.l.A.A.  Journal,  v.  2,  Oct.  1964,  p.  1788-1794. 

9.  TSUI,  TIEN-YU,  and  LAKSHM1KANTHAM,  C.  The  Nonlinear  Dynamic  Response  of  Stiffened  Shells  Under 
Compressive  Axial  Impulse.  Presented  at  the  4th  Canadian  Congress  of  Applied  Mechanics,  Montreal,  Canada, 
1973. 

10.  LAKSHMIKANTHAM,  C.,  and  TSUI,  TIEN-YU.  Dynamic  Stability  of  Axially  Stiffened  Imperfect  Cylindrical 
Shells  Under  Axial  Impulsive  Loading.  Accepted  Publication  by  A1AA  Journal. 
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The  assumptions  basic  to  the  linear  problem  are  similar  to  those  of 
Reference  6:  a)  the  stiffeners  are  close  enough  that  they  can  be  ''smeared" 
over  the  shell;  b)  the  stiffeners  are  essentially  beam-like  elements 
contributing  to  the  membrane  and  bending  resistance  of  the  shell  but 
not  to  the  shear  resistance;  c)  the  entire  shell  including  the  stiffeners 
is  activated  in  the  buckling;  that  is,  local  instability  failures  are 
ruled  out. 

Figure  1  describes  the  basic  geometric  parameters  of  the  shell. 

With  the  normal  to  the  shell  surface  taken  as  positive  inward,  positive 
(negative)  values  of  Es,  Ef  indicate  inside  (outside)  location  of  the 
stiffeners  with  respect  to  the  shell. 

The  unstiffened  shell  is  characterized  by  three  parameters  L,  R,  3' 

H  (length,  radius,  and  the  thickness  of  the  shell)  which  can  be  combined 
into  the  so-called  Batdorf  parameter  Z  =  L2/RH,  the  ranges  of  whose 
values  are  used  to  describe  a  "short"  or  a  "moderate  length"  cylinder. 

The  effective  increases  in  the  area  of  section,  the  section  area 


Figure  1.  Shell  Geometry 
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moment,  and  the  section  torsional  stiffness  of  the  shell  due  to  the 
presence  of  stiffeners  are  best  characterized  by  nondimen sional  parameters 
a,  n  and  y  with  subscripts  s  or  r  to  denote  a  stringer  (axial  stiffener) 
or  a  ring  (frame  or  circumferential  stiffener).  These  parameters  are 
the  ratios  of  the  properties  of  the  stiffener  section  relative  to  a 
unit  shell  section  between  the  stiffener  center  lines.  We  have  assumed 
that  the  stiffeners  have  the  same  material  properties  as  those  of  the 
shell;  otherwise  the  relative  stiffener-to-shell  moduli  ratios  are  to 
be  included  in  the  definition  of  a,  n  and  y. 

Although  these  parameters  are  not  independent  of  each  other,  it  is 
still  convenient  to  consider  them  so,  from  a  computational  point  of 
view.  For  slender  and  deep  stiffeners  of  rectangular  sections  (such  as 
is  assumed  here)  one  can  express  y  in  terms  of  a  multiplied  by  a  constant. 
Hence  the  chief  stiffener  parameters  of  interest  are  a,  n  together  with 
E  the  eccentricity  of  the  stiffener. 

In  order  to  derive  the  basic  equations,  we  introduce  the  following 
nondimensional  quantities: 

All  lengths  (including  the  coordinate  distances  X,  Y)  are  normalized 
with  respect  to  R,  the  shell  radius.  Thus  we  write: 

[u,  v,  w,  x,  y,  eg,  er]  =  [U,  V,  W,  X,  Y,  Eg,  E^J/R  (1) 

Dimensionless  stress  and  moment  resultants  are  given  by: 


n  ,  n  ,  n  ,  n  =  [N  ,  N  ,  N  ,  N  ]/B  B  =  Eh/(l-v2) 

x’  y’  xy  x  1  x’  y*  xy’  xJ 


m.m.m  ,  m  =[M,M,M  ,M  1/D  D  =  Eh3/12(l-v2) 
x’  y’  xy’  yx  1  x’  y’  xy’  yxJ  1 


(2) 
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The  linear  strain-displacement  relationships  at  the  middle  surface  of 
shell  are: 


c  =  u ,  ;  e  =  v ,  -  w  ;  e  =  1/2  (u,  +  v>,J 

x  *x  y  y  xy  y  x 

Rtc  =  w,  ;  Rtc  =  w,  •  Rtc  =  w,  =  Rtc 
x  xx  y  ’yy  *  xy  xy  \  x 


Tlie  stress-resultant  strain  relationships  are  given  by: 


(3) 

(4) 


The  nonzero  submatrices  M2  and  in  Eq.  (5)  provide  the  coupling 
effects  of  the  stiffener  eccentricity. 

In  Eq.  5  Ss  (8r)  is  not  a  new  parameter  but  is  related  to  ag  (a^) 
through  D3g  =  Bas> 


The  governing  equations  are  then  given  by: 


n  +n  =  n  +n  =0 
x,x  xy,y  xy,x  y,y 


(6a,  6b) 


1  H  L_  [wj  +  a  e  e°  +  a  e  c°  -n  w,  +  n  =  0  (6c) 

y-,  r2  3  1  J  s  s  x,  xx  r  r  y,yy  x  'xx  y 


s  s  x,  xx  r  r  y,yy  x  xx 
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where  is  a  linear  differential  operator  given  by 

L,  (  )  =  (1  +  n  )  (  ),  +  2  {1  +  (Y  +  Y  )}  (  ),  + 

3  v  s'  xxxx  2  v  s  r  v  "xxyy 


0  *  V  (  5 'yyyy 

and  nx  is  the  (dimensionless)  external  compressive  load. 


(6d) 


Eq.  (6c)  reduces  to  the  standard  (linear)  form  when  all  stiffness 

and  eccentricities  are  neglected. 

Eqs.  (6a,  6b)  are  identically  satisfied  by  introducing  a  stress 

function  of  such  that  n  =  f,  :  n  =  f,  and  n  =  -f. 

x  yy  y  xx  xy  xy. 

Making  use  of  the  stress-strain  relationships  Eq.  (5)  and  introducing 

the  stress  function  f,  Eq.  (6c)  may  be  rewritten  as: 

"17  F  L3  (w)  ~r-  L2  (f)  +  F~  L4  (w;)  "  "x  Wjxx  +  f’xx  "  0  (7) 

rs  rs 


where  L 2  and  L 4  are  linear  operators  given  by 

L~  (  )  =  va  e  (  ),  -{(]  +  a  )  a  e  +(l+a)aF)(  ), 

2  s  s  v  ’xxxx  v  s  r  r  v  s  s 

va  e  (  ) , 

r  r  yyyy 


xx/y  + 
(8a) 


and  L.(  )=(l+a)(ar)  (  ),  -2vaa  e  e  (  ), 

4  r  3  s  v  ’xxxx  s  r  s  r  J,xxyy  + 


(1  +  a  )  (a  e  )  (  ), 


r  r 


yyyy 


(8b) 


and  A  =(l+a)(l+a)-  v2 
rs  s  r 

A  second  independent  equation  involving  f  and  w  is  provided  by  the 

compatibility  reiat’  xn  which  is  easily  deduced  from  Eq.  (3)  to  be 

c°  +  e°  -  2e°  =-w,  (Q) 

x,yy  y,xx  xy,  xy  xx  J 

Eq.  (91  can  be  written  in  terms  of  f  and  w  by  again  making  use  of  Eq.  (5). 

This  results  in 
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(10) 


L.  (f)  =  L.  (w)  -  A  w, 

1  v  '  2  rs  xx 


rs 


where  L  (  )  =  (1  *  a  )  (  ),  +  2  (-j - v  )  (  ),  +  (!  +  “  ) 

i  s  xxxx  i  -  v  xxyy  r 

(  }’yyyy  (ii) 

Thus  the  basic  equations  of  the  static  buckling  problem  are  Eqs. 
(7)  and  (10). 


Now  we  notice  that  the  operators  L^,  L^,  L^,  are  linear  with 
constant,  coefficients  (shifted  biharmonic  types). 

>ince  the  buckling  modes  are  essentially  harmonic,  w  and  consequently 
f  are  both  harmonic.  Hence  we  see  that  the  following  operations  with 
the  operators  hold: 


L.  [L.  (♦)]  =  Lj  [L.  (♦)]  (12) 

and  also  [  (  ) >xx]  =  [Li  (  )]>xx  i,  j  -  1,  2,  3,  4  (13) 

Thus  we  can  eliminate  f  from  Eq.  (7)  by  operating  with  upon  Eq.  (7) 
and  as  a  result  have 


12 


2 

£  2L1  4  <"> 


r~  4  4  <«  *  r~  4  4  <“>  -  4  4  <“-xx>  * 


L,  f. 

1  xx 


rs 
=  0 


(14) 


By  making  use  of  Eqs.  (12),  (13)  and  (10)  we  finally  have 
2 

12  Ip  4  4  -  r-  h  l4  <“>  -  4  4  <">]  4  A™  "-xxxx  -  2  4<"-xx> 


rs 


*  \  q^-xx5 


o 


(15) 


Eq.  (15)  is  thus  the  Donnell  equation  for  the  eccentrically  stiffened 
shell  problem  with  only  w  appearing  in  it. 
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We  seek  a  solution  of  the  above  by  choosing  the  following  form  for 


w: 

w  =  w  cos  mirx  cosny 

where  it  =  -nR/L 

Using  Eq.  (16)  in  Eq.  (15)  we  find 


(16) 


n  R 

JL  _  in  = 

A  H 
rs 


rs 


(mir)  1  C 


-  r  rs  - 


rs 


!  ,  D 

_  {  rs  - 

2 

B  \ 
rs 

P 

z 

1 

(17) 

2  8 

+ 

_  + 

rs 

A 

rvl 

*= 

B 

^ _ / 

P 

A  P 

rs 

rs 

rs 

rs  rs 

(i*  ('-2v> 

(Ys  * 

y  ))  e2+ 

r 

(l  +  nr) 

84 

2  2 
D  =  (a  e  1  (1  +  a  )  -  2v  aaee  82  +  (1  +  a  )  (a  e  )  B4 
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The  buckling  load  is  obtained  by  minimizing  Eq.  (17)  with  respect 
to  the  wave  numbers  n,  m. 

In  deriving  Eq.  (17),  we  notice  that  no  emphasis  has  been  placed 
on  the  boundary  conditions.  This  is  as  it  should  be.  It  is  inherent 
in  the  classical  shell  stability  heory  that  the  eigenvalues  are 
obtained  without  reference  to  boundary  conditions. 

DISCUSSION  OF  THE  BUCKLING  EQUATION 


Eq.  (17)  constitutes  the  fundamental  equation  from  which  the  (minimum) 
buckling  load  is  obtained  for  a  given  shell  with  prescribed  geometry  and, 
furthermore,  it  can  be  demonstrated  that  it  is  identical  with  the  equation 
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used  by  Singer  et  al.6  We  now  proceed  to  deduce  several  interesting 
results  based  on  the  form  of  Eq.  (17)  together  with  the  stiffening  para¬ 
meters  defined  in  Eq.  (18). 

Firstly,  we  notice  that  the  eccentricity  parameters  es,  er  appear 

explicitly  only  in  the  terms  B  and  D  .  Thus  B  and  D  are  the 
v  1  1  rs  rs  rs  rs 

coupling  effects  due  to  eccentricity.  If  the  eccentricities  were  ignored, 
(either  owing  to  symmetrical  stiffeners  or  as  a  higher  order  effect) 


Brs  and  Drs  would  drop  out,  and  the  resulting  expression 


1  (imr) 


rs  + 


12  Z 


would  admit  positive  (and  hence  physically  admissible)  solutions  for  all 
choices  of  a  and  n. 

However,  B  involves  first  powers  of  e  ,  e  while  D  is  quadratic 
in  them.  Hence  the  sign  of  the  last  term  in  Eq.  (17)  is  dictated  by  the 
sign  of  the  eccentricity  terms.  This  brings  out  a  distinctive  feature 
of  the  problem  of  an  eccentrically  stiffened  shell.  The  presence  of  a 
minus  sign  within  the  square  brackets,  together  with  the  possibility  that 
the  last  term  could  be  negative,  makes  only  certain  combinations  of  the 
stiffening  parameters  yield  physically  admissible  solutions.  It  is  not 
hard  to  see,  given  the  structure  of  Eq.  (17),  that  there  are  admissible 
combinations  of  the  parameters  yielding  negative  or  zero  values  for  the 
buckling  load!  It  is  even  more  remarkable  that  the  same  set  of  physical 
parameters  which  would  yield  admissible  (positive)  solutions  based  on 
Eq.  (19)  would  cease  to  do  so  if  the  symmetry  of  the  stiffeners  is 
distorted. 
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A  second  observation  which  we  can  make,  and  to  which  we  shall  come 
back,  is  the  absence  of  Z  explicitly  in  the  last  term  of  Eq.  (17),  which 
makes  certain  solutions  valid  fcr  all  Z  ranges. 

Eq.  (17),  involving  both  rings  and  stringers,  is  applicable  to  a 
general  (grid)  stiffened  shell.  Particular  cases  of  axial  or  ring 
stiffeners  are  obtaii.^d  from  Eq.  (17)  br  setting  the  appropriate  stif¬ 
fening  terms  equal  to  zero.  A  grid  stiffened  shell  is  essentially  a 
composite  of  axial  and  ring  stiffened  cylinders,  and  hence  its  behavior 
is  understood  best  by  studying  the  axial  and  the  ring  cases  separately. 
This  paper  confines  its  attention  to  the  axial  and  ring  cases  only. 

AXIALLY  STIFFENED  SHELL 


Ail  the  parameters  with  subscript  r  are  set  equal  to  zero  in 
Eqs.  (17)  and  (18)  and  after  some  manipulation  we  can  obtain  the 
governing  equation  for  the  axially  stiffened  shell  under  axial  compres¬ 
sion  as  follows; 
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a  constant  w  displacement  which  is  of  no  interest  to  the  buckling  problem). 
The  lowest  of  these  loads  for  all  m,n  combinations  is  the  buckling  load 
for  the  given  shell. 

From  a  theoretical  point  of  view,  however,  it  is  more  convenient 
to  assume  ab  initio  a  buckling  mode  and  investigate  the  range  of  its 
validity.  We  first  consider  the  axisymmetric  case  (n=0) . 


Axl symmetric  Case 

This  is  obtained  by  setting  6=0  in  Eqs.  (.TO)  and  (21).  The 
resulting  equation  is  then  minimized  with  respect  to  m  to  yield  the 
following : 


r  *  f  1  *  "s  _  (vs)2n 
(r  rrr^  L-rr-  tt^I) 


/  2va  e 
l/2  s  s 

+  - 7-r- 


A  (1  +  a  ) 
s  s 


The  sign  cf  the  last  term  in  Eq.  (22)  denends  upon  the  sign  of  eg. 
We  can  thus  distinguish  the  two  cases  of  inside  ( +  ve  e^  an^  outside 
(-  ve  es)  stiffeners. 


We  notice  that  for  real  values  of  n  . , 

axi 

ns  -  12  (egas)?  /  (1  ♦  ag)  -  1  s  n  *  (23) 

should  hold  in  Eq.  (22).  It  is  interesting  to  note  that  for  ng  greater 
★ 

than  n+  such  that  the  square  root  term  is  positive,  the  predicted  load 
is  affected  by  the  sign  of  the  last  tern.  Thus,  for  e^  positive  the 
load  is  higher  than  that  with  eg  negative.  Hence  the  inside  stiffener 
case  would  yield  higher  axisymmetric  loads  than  outside  stiffeners. 

It  follows  that  for  negative  es>  the  minimum  load  is  zero  and  the  corre¬ 
sponding  n  value  is  given  by: 
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■  v^r.  t  •** . »?*%*» "  ’■ony^TOw  »r^|y.i!w.c,-iwA-  j,  •«*»«.• 
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ns  «  n_  *  12  (es«s)2  /  Ag  -  1  Jdth  nflxi  (-)  =  0)  (24) 

*  * 

n  is  greater  than n+ as  Ag  (il  +  as  -vz)  <  1  +  ag  for  positive 

v.  However  the  difference  between  the  two  n  values  is  very  little. 

For  v°ry  large  n  values  the  difference  is  of  the  order  e^2. 

An  important  feature  of  the  axisymmetric  solution  is  that  it  is 
independent  of  Z.  Hence  if  the  solution  is  valid  it  will  hold  for  both 
'short'  ana  'moderate-length'  cylinders. 

Asymmetric  Case 

The  asymmetric  solution  is  obtained  from  Eq.  (20)  as  the  lowest 

possible  load  for  all  possible  integer  combinations  of  m  and  n  (n>0) 

for  all  given  values  of  a,  e,  n,  and  Z.  Obviously,  a  large  number  of 

computations  has  to  be  performed  to  cover  a  wide  range  of  these  parameters. 

However,  a  slightly  modified  form  of  Eq.  (20)  is  more  convenient  for 

theoretical  comparison  with  the  axisymmetric  solution  Eq.  (22).  This 

approximate  form  is  obtained  by  noting  that  for  the  axisymmetric  case  of 

a  moderate  length  cylinder,  generally  n>>m  and  hence  62  can  be  taken 

much  larger  than  unity  in  Eq.  (20)  and  thus  we  write: 

1  (mir) 2  C  2  (a  e  )2  82  Z  1  2ae 

n  =_  _  l  _  -  _  ]  +  _  _  - _ 

A  Z  12  P  (mu)2  P  A  P  (25) 

s  s  sss 

Even  further  simplification  is  possible  if  we  compare  the  sizes  of 

4 

some  of  the  terms  with  8  in  Eqs.  (20)  and  (21).  Generally,  a  and  e  have 
numerically  the  smallest  range  of  values,  seldom  exceeding  a  number 
such  as  10  for  practical  cylinders;  however,  n  has  a  wide  range  and 
can  have  values  of  the  order  104.  Thus  retaining  terms  of  the  order 
B4  in  Eqs.  (20)  nd  (21)  we  can  rewrite  Eq.  (25)  as: 
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.  1  (mir) 2  l  +  n+S4  Z  1  ?  a  e 

n  = _ [  5  ]  * _ -  5  5  (26) 

A  Z  12  (mir)2  64  A  &2 

s  '  s 

Certain  features  of  the  axi symmetric  solution  in  comparison  with 
the  axisymmetric  solution  Eq.  (22)  became  quite  apparent  from  the 
form  of  Eq.  (26).  Firstly,  we  notice  that  the  iign  of  the  last  term 
in  Eq.  (26)  is  opposite  to  that  in  Eq.  (22).  Thus  the  stiffener  effects 
is  reversed  in  the  asymmetric  case:  buckling  load  for  es  positive 
(inside  stiffeners)  is  lower  than  that  for  negative  es>  Secondly, 
the  solution  of  Eq.  (26)  is  strongly  Z-dependent.  Also  in  Eq.  (26) 
there  is  no  apparent  lower  limit  to  n  values  in  contrast  to  the  n* 
values  of  tqs.  (231  and  (24);  hence  in  the  regions  where  the  axisym¬ 
metric  solution  is  inadmissible,  the  asymmetric  solution  will  prevail. 
Consequently,  some  valid  comparisons  can  be  readily  made  between  the 
two  modes;  however,  it  is  convenient  to  distinguish  the  two  cases  of 
eccentricity. 


Outside  Stiffeners  (es  negative)  :  Comparison  of  Solutions 

For  negative  es>  all  the  terms  of  Eq.  (26)  are  positive,  whereas 
the  minimum  n  .  is  zero  from  Fq.  (24).  Hence  we  suspect  that  for 

HX 1 

*  > 

any  a,e  combination,  in  the  vicinity  of  n  (i.e.  ns  -  n*)»  given 
from  Eq.  (24),  the  axisymmetric  solution  will  prevail.  For  values  of 


ng  much  larger  than  n*,  Eq.  (22)  shows  a  monotone  behavior  and  for 

very  large  ns  values,  Eq.  (22)  will  approach  the  limit. 

1  +  n  i/2  1  1/2 

n  .  =  Z  ( - )  ( - 1 

axi  .  '  ,  _  , 

1  +  a  12  A 

s  s 
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However  for  large  ns  values,  Eq.  (26)  is  still  modified  by  Z  and 
other  factors.  Intuitively,  we  expect,  then,  that  the  asymmetric  solu¬ 
tions  will  prevail  for  large  values.  Again  for  ns  lower  than  the 
critical  value  for  axisymmetric  solution  (=n*) ,  the  asymmetric  solution 
being  admissible,  we  can  assume  that  we  will  have  the  asymmetric  solution. 
Thus,  in  an  ns  -  as  plane  (fcr  fixed  values  of  es),  we  expect  two  zones 
of  asymmetric  modes  interspersed  by  a  (narrow)  zone  where  axisymmetric 
solution  prevails. 


Inside  Stiffeners  (e_,  positive)  :  Comparison  of  Solutions 

Here,  from  Eq.  (26)  we  see  that  n  values  can  decrease  to  zero 

★ 

whereas  the  lowest  possible  axisymmetric  load  (occurring  at  =  n+) 
is  seen  from  Eq.  (22)  to  be 


n  .  (+)  = 
-axi  v  1 


zva  e 
s  s 


(28) 


A  (1  +  a  ) 
s  v  s 


Hence  we  can  expect  that  there  may  be  no  axisymmetric  zone  at  all 
for  the  Z  values  in  the  moderate  length  region. 


Numerical  Results:  Discussion  of  Zones 

In  order  to  give  substance  to  the  above  arguments,  a  number  of 
a,  e,  n  combinations  were  used  to  study  Eqs.  (20)  and  over  a  wide  range 
of  Z  values  from  102tol04  a  was  varied  from  0  to  10;  e  was  taken  as 
+  1,  3,  5,  7,  and  10.  For  each  combination  of  the  parameters,  a  minimum 
of  Eq.  (20)  was  found  by  considering  integer  values  of  m  (1  to  20)  and 
n  (1  to  50)  on  a  Univac  1106  computer. 


I 

3 
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The  results  are  shown  in  Figure  2.  For  every  cace  the  n-a  diagram 
shows  four  typical  zones:  Zone  I  of  asymmetric  m=l  mode;  Zone  II  of  Z- 
dependent  axisymmetric  and  asymmetric  (m=l)  modes:  Zone  III  of  axisymmetric 
mode  and  Zone  IV  of  asymmetric  mode  where  generally  m  ^  1. 

Zone  I  occurs  for  Z  generally  in  the  moderate  length  range,  though 
this  range  decreases  with  increasing  e.  Thus,  for  the  highest  e 
considered,  i.e.  e  =  10,  the  Z  range  of  Zone  I  was  from  150  to  104  while 
for  e=l,  Z  could  be  as  low  as  20. 


0.1  0  2  0  4  1.0  2.0  4.0  10.0  0.1  0.2  0.4  1.0  2.C  4.0  10.0 
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Zone  II  is  a  transition  zone  where  axi symmetric  solutions  go*em  for 

Z  values  which  increase  with  decreasing  n,  until  Zone  III  values  of  n 

are  reached,  where  axisymmetric  solutions  givems  all  Z  values.  Generally 

Zones  II  and  III  form  a  narrow  band  and  the  curve  distiguishing  the 

two  becomes  asymptotic  to  the  curves  of  Zone  I  and  Zine  III.  The  lower 

bound  of  Zone  III  is  the  n*  curve  defined  in  F.q.  (24). 

For  positive  e  values  also, there  exist  the  same  four  zones  though 

the  middle  zones  are  even  less  prominent  than  in  the  negative  e  case. 

It  is  interesting  to  note  that  the  same  curve  forms  the  lower  bound  of 

Zone  I  in  both  the  positive  and  negative  cases.  It  may  also  be 

mentioned  that  the  values  of  shell  parameter  studied  by  Singer  et  al6 

fall  within  Zone  I  for  the  axially  stiffened  case. 

Summarizing  the  axially  stiffened  case,  we  have  for  every  a,e 

combination  a  minimum  n  above  which  the  buckling  is  in  asymmetric  m=l 

mode,  valid  for  t'.e  shell  parameter  Z  with  a  minimum  range  of  150  to  104; 

for  lower  n  values,  axisymmetric  solutions  prevail  up  to  a  critical 
★ 

value  of  n  below  which  once  again  asymmetric  modes  (though  m  ^  1  now) 
occur . 

From  a  predictive  point  of  view,  where  we  wish  to  design  a 

cylinder  to  buckle  in  the  m=l  mode,  we  are  obviously  interested  in  the 

bounding  curve  of  Zone  I.  The  following  approximate  expression  was 

found  to  provide  a  close  upper  bound  on  the  curve: 

n  =  18  Cases)2  (29) 

1  +  a 

s 

valid  for  a'.  1  the  ranges  of  the  parameters  tested. 
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RING  STIFFENED  SHELL 


For  the  ring  stiffened  shell,  we  set  all  terms  with  subscript  s 
equal  to  zero  in  Eqs.  (17)  and  (18)  and  after  some  manipulation  obtain 


n  = 


(mu)2 

A  Z 
r 


[*■ 


(a  e  ) 
r  r 


2  S4  ( 


+  262  i 

:->] 


(mir) 1  Pv 


2a  eg2  (l-v82) 
r  r 

A  P 
r  r 

(30) 


where 


Cr  =  1  ♦  2  (1  ♦  yr)62  ♦  (1  *  nr)  84 

Pr  =  1  +  2  (1  ♦  ar/l-v)82  +  (1  ♦  ar)B4  (31) 

A  =  1  +  a  -v2  8  =  n/mir 

r  r 


Axisyncnetric  Solution 

The  axisymmetric  solution  is  readily  obtained  by  setting  8=0  in  Eqs. 
(30)  and  (31)  and  minimizing  the  resulting  expression  with  respect  to  m. 
Thus  we  have: 

nax.  =  2/(12  Ar)1/2  i  1/ [3 (1  ♦  ar  -  v2)]1/2  (32) 

Notice  that  this  solution  is  independent  of  both  n  and  e. 


Asymmetric  Case:  A  Closed  Form  Solution 

The  asymmetric  case  is  based  on  Eq.  (30)  by  the  usual  minimization 
with  respect  to  m,  n  (n  >  0) .  However,  it  is  possible  to  obtain  a 
simple  closed  form  solution  for  the  particular  case  of  m=l  and  62  >  >  1. 
Thus  by  retaining  terms  of  che  order  B4  in  Eqs.  (30)  and  (31),  we  get 
the  following  approximate  expression: 
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1  ♦  nr  (a  e  )2 

_  q4  r _ £  r  r 

A_Z  P  1  12  1  +  a 


1  + 


2va  e 
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(1  +  a  )8'1  A  (1  +  a  ) 
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(33) 


which  yields  the  following  minimum 


£rV\,l/2  ^Vr 

1  T  “r  +  Ar  (1  +  ar) 

(34) 

We  notice  that  the  form  of  Eq. 
(34)  is  identical  with  the  axisymmetric  solution  for  the  axially 
stiffened  case,  Eq.  (22).  This  solution  again  leads  to  possible  zero 
value  for  negative  e^  (outside  stiffeners).  Hence  we  can  suspect  that 
for  nr  values  very  close  to  (and  greater  than) 


1  1 


1  +  n. 


*■2(1- 


A  (1  -  a  )  1  12 

r  ‘  r 


It  is  instructive  to  examine  Eq.  (34) 


r  (-) 


12  (a  e  ): 
-  r  rJ 


-  1 


(35) 


r 

the  asymmetric  m=l  mode  will  prevail  for  at  least  the  negative  er  case. 

★ 

Furthermore,  for  nr  <  ir  ^  ^  we  see  that  Eq.  (34)  will  not  hold  as 

*  * 

it  will  yield  negative  values  for  n;  hence  for  nr  <  nr  we  should  expect 
other  modes,  more  general  asymmetric  or  the  axisymmetric  modes,  to 
prevail . 

Before  we  try  to  compare  the  solutions,  we  can  examine  the  conse-* 
quences  of  letting  m  be  of  the  same  order  as  n  in  Eq.  (30);  the  essential 

result  is  that  82  <  1  and  the  last  term  in  Eq.  (30)  gets  modified  to 
„  (m*)2  r  Cr  84  +  2B2  Z  1  2arerB2 

n  =  "7TZ  [  12  '  (arV2  g4  {  Y-v  }  J+  liinfy2  F"  '  (36) 


P 

r 

which  is  additive  for  negative  e^  and  subtractive  for  positive  er< 
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Comparison  of  Solutions:  Outside  Stiffeners  (negative  e  ) 

Of  the  three  possible  solutions,  we  can  anticipate  that  for  valuer 

* 

of  n  slightly  in  excess  of  nr  given  by  Eq.  (35),  Eq.  (34)  will  provide 

★ 

the  lowest  possible  values;  when  n  is  substantially  larger  than  n  , 

the  monotonic  behavior  of  Eq.  (34)  suggests  that  the  only  possible 

solutions  will  be  the  axisymmetric  case  of  Eq.  (32)  or  the  more  general 

asymmetric  case  corresponding  to  Eq.  (36).  However,  here  again  since 

the  last  term  is  additive  for  negative  e  ,  and  also  that  Cr  will  be 

quite  dominant,  one  can  intuitively  opt  for  the  axisymmetric  case  which 

* 

is  independent  of  n  and  Z.  For  values  of  n  below  n  ,  the  only  solutions 
in  contest  are  the  axisymmetric  or  the  general  asymmetric.  Here  again, 
since  Cr  may  itself  be  small,  the  terms  in  the  square  bracket  may  be 
substantially  low  to  warrant  the  prevalence  of  the  general  asymmetric 
solution  over  the  axisymmetric  solution. 

Comparison  of  Solutions:  Inside  Stiffeners 

For  inside  stiffeners,  with  er  positive,  Eq.  (36)  could  give  us 
substantially  lower  values  as  the  last  term  is  negative,  compared  to 
the  m=l  solution,  Eq.  (34),  where  the  last  term  is  positive.  Hence 
the  probable  modes  are  the  axisymmetric  and  the  general  asymmetric 
modes.  A  further  careful  examination  of  Eq.  (36)  shows  that  it  is 
strongly  dependent  on  Z,  generally  increasing  with  Z  (as  £2  <  1  and 
all  the  terms,  such  as  C  ,  tend  to  a  limit).  Hence  for  high  Z 
values  Eq.  (36)  may  yield  high  values  of  n  which  are  higher  than  the 
axisymmetric.  For  lower  values  of  Z,  Z/ (mu) 2  will  be  substantially 
smaller  and  the  sizes  of  the  positive  and  negative  terms  will  be  of  the 
same  order;  thus  yielding  quite  low  values. 
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Numerical  Results  and  Discussion 


Eq.  (30)  was  studied  for  several  combinations  of  a,  e,  and  n 
over  a  range  of  Z  from  102  to  104  and  m,n  values  in  the  range  (1  to  50) 
each.  Eq.  (32),  being  a  function  of  a  only,  was  computed  for  all  the 
a  values  of  interest. 

The  role  on  stiffener  eccentricity  on  Eq.  (30)  was  found  to  be 
that,  for  negative  e  (outside  stiffeners),  n  (i.e.,  the  asymmetric  load) 
was  relatively  independent  of  Z,  depending  only  on  n for  a  shell  of 
given  a,  e. 

However,  for  positive  e,  n  was  strongly  dependent  on  both  Z  and  n. 
Since  this  difference  in  behavior  distorts  the  regions  of  asymmetric  and 
axi symmetric  behavior  it  is  convenient  to  treat  the  two  cases  of 
eccentricity  separately. 

a)  Outside  Stiffeners:  The  solutions  of  Eq.  (30)  for  negative  e  were 
almost  independent  of  Z  (in  the  range  or  interest)  and  increased  rapidly 
with  n  so  that  for  every  a,  e  case  it  was  easy  to  find  a  minimum  n 
above  which  tbs  asymmetric  solution  of  Eq.  (30)  was  always  higher  than 
the  axisymmetric  solution. 

Figure  (3)  shows  a  typical  n-a  diagram  for  a  suitable  |e|  value. 

The  diagram  consists  of  twi  major  zones  of  axisymmetric  and  general 

asymmetric  modes.  The  dividing  curve  happens  to  coincide  with  the 
★ 

rir  (  j  curve  given  in  Eq.  (35),  which  is  based  on  an  m=l  solution. 

★ 

This  mode  is  essentially  confined  to  the  n  curve;  for  values  of  n 

★ 

even  slightly  larger  than  n  we  get  only  axisymmetric  solution  as  the 
buckling  mode. 
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Hence  from  a  predictive  point  of  view  any  value  of  n  slightly 
* 

larger  than  nr  ^  ^  would  give  us  the  minimum  combinations  for  an 
axisymmetric  buckling. 

b)  Inside  Stiffeners:  For  positive  e,  the  dependence  of  Eq.  (30)  on  Z 
exhibits  itself  in  a  remarkable  way.  For  Z  <  103  generally  n  was 
independent  of  Z,  but  in  the  range  10 3  >  Z  >  104  it  almost  increases 
monotonically  with  Z.  This  behavior  is  generally  true  for  low  n  as 
well  as  high  n  as  shown  in  an  example  in  Figure  4.  (for  a  combination  of 
a,  e  which  could  be  taken  as  light  stiffening).  The  axisymmetric  line  (being 
independent  of  Z,  n>  e)  is  also  sbjwn  for  the  a  considered.  Thus  it 
becomes  clear  that  for  Z  >  103  we  can  have  axisymmetric  buckling 
(provided  we  have  a  minimum  n  required) .  For  every  case  studied  we 
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Figure  4.  Variation  of  Buck*  >  3  Loads  Versus  Z  for  Cylinder  with  Inside  Radial  Stiffeners 

(a,  *  0.5,  er  =  3) 

found  that  for  Z  >  103,  we  could  find  acceptable  n  values  beyond 
which  axi symmetric  solutions  become  possible.  It  is  also  clear  from 
Figure  4  that  the  higher  the  Z  value  the  minimum  n  for  axisymmetric 
buckling  becomes  lower.  Figure  5  shows  the  minimum  n  values  for  a 
shell  of  Z  =  2000. 

The  dividing  curve  between  the  axisymmetric  and  the  asymmetric 
solution  in  Figure  5  can  be  represented  by 


n  =  C 


12  (a  e  )2 
1  +  a_ 


(37) 


where  the  values  of  C  are  given  in  Figure  6  for  different  values  of  e. 
Returning  to  Figure  4,  we  notice  that  for  Z  <  103  the  axisymmetric 
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Figure  5.  Governing  Buckl.ng  Modes  for  Radially  Stiffened  Cylinder  Under  Axial  Compression 


solution  was  independent  of  Z,  but  that  for  the  particular  n  shown  (200) 
it  was  slightly  below  the  axisymmetric  value.  However,  no  matter  how 
high  the  value  of  n  was  chosen,  the  asymmetric  solution  for  Z  <  103 

did  not  exceed  the  axisymmetric  line.  Thus  a  second  remarkable  feature 
of  the  positive  eccentricity  emerges:  for  Z  <  103  the  axisymmetric 

solution  seems  be  almost  impossible.  The  reason  is  probably  that 
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Figure  6.  Graphs  for  Constant  C  in  Equation  (37) 
for  Z  <_  103  the  buckling  is  a  'many-bayed'  phenomenon  with  each  bay 

buckling  as  a  short  cylinder.6  However,  a  point  remains:  the  axisymmetric 

solution  is  certainly  possible  for  Z  >  10 3  which  Reference  6  seems  to  regard 

as  impossible  for  any  Z. 

A  further  remark  is  also  relevant:  the  axisymmetric  solution, 
being  so  close  to  the  asymmetric  solution  (at  least  for  the  q  shewn  in 
Figure  4),  provides  a  close  approximation  for  the  buckling  load  for  102 
<  Z  <  103  for  a  great  number  of  a,  e,  q  combinations. 

Summarizing  the  ring-stiffened  case,  both  axisymmetric  and  asymmetric 
buckling  are  possible  and  for  outside  stiffeners  it  occurs  for  all  Z  ranges 
while  for  inside  stiffeners  it  is  confined  to  moderate-to-long  Z  ranges. 
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CONCLUSION 


The  primary  aim  of  this  paper  has  been  to  concentrate  on  the  effect 
of  stiffener  eccentricity  on  the  buckling  modes  rather  than  on  the 
buckling  loads  ;  to  present  as  simply  as  possible  the  parameter  combi¬ 
nations  for  which  a  particular  mode  would  prevail.  In  some  cases  we 
have  succeeded  in  presenting  a  simplified  relationship  between  the 
parameters  which  would  give  minimum  parameter  combinations  for  moderate- 
length  (102  <  Z  <  104)  cylinders  buckling  in  the  dominant  mode.  By 

concentrating  on  the  buckling  equations  (within  the  assumptions  of  the 
theory)  we  have  avoided  some  of  the  conclusions  made  in  References  4  or  6, 
i.e.,  "that  axisymmetric  buckling  occurs  for  'short'  cylinders'  or  "that 
axisymmetric  buckling  cannot  occur  for  ring-stiffened  shells  with  inside 
stiffeners". 

This  a  priori  knowledge  of  buckling  modes  is  extremely  useful  in 
the  computationally  efficient  integration  schemes  for  the  related  problem 
of  the  dynamic  stability  of  stiffened  shells. 
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